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Introduction 

The physics of the scalar mesons has been challenging 
for decades. Scalar mesons are difficult to detect, since 
they decay isotropically and have, in general, large widths. 
There are many candidates with mass bellow 2 GeV/c^. 
Some states are now well established, while others remain 
controversial. In any case, there is a large overlap between 
states in this region of the spectrum. 

An additional problem arises from the richness of the 
low energy strong dynamics, allowing other 0++ configu- 
rations than the usual qq. Although none of these 'exotic' 
configurations has been clearly established, some would 
populate the 1-2 GeV/c^ region, mixing with the regular 
qq mesons. The identification of the qq scalar nonet(s) is, 
thus, a rather complicated task, which can only be accom- 
plished if one combines data from different types. 

This paper is focused on two issues: the Ktt spectrum 
near threshold - the kappa problem -, and the vrvr spectrum 
between 1.2 and 1.5 GeV/c^ - the /o(1370) problem. Both 
problems will be analysed from the perspective of the D 
meson decays, with some additional information from t 
and B decays. 

Hadronic and semileptonic decays of D mesons have 
unique features that make them a key to light meson spec- 
troscopy, in particular to the study of the Ktt and tttt scat- 
tering amplitudes in S-wave. First, these are the only pro- 
cess that allows us to access the Ktt/tttt spectrum continu- 
ously from threshold up to 1.5 GeV/c^. In Kn and tttt 
scattering, the production of scalar resonances near thresh- 
old is suppressed by the Adler zeroes. No such effect is 
observed in D decays. In scattering, the large nonreso- 
nant component forms a continuum background on top of 
which the scalar resonances are found. The interference 
between the broad states and this continuum distorts the 
resonance line shape, and is always a difficult problem to 
be accounted for In D decays, the nonresonant compo- 
nent is usually small. In D decays, channels with two 
identical pions in the final state have a largely dominant 
S-wave component. Decay modes like — > K^tt^tt^ 
and _D+, 13+ 7r^7r+7r+ are easy to be reconstructed and 
have large branching fractions. There are plenty of good 
data with very low background on these 'golden' modes. 
In the first year of the LHCb operation, we will enter in the 
regime of 'infinite' statistics. 

The decay of a Z3 meson is, obviously, a very complex 
process. It is initiated by the c ^ s{d) weak transition. 
This transition is embedded in a strongly interacting sys- 
tem, in the non-perturbative regime, from which the final 
state hadrons emerge. No precise quantitative description 
based on first principles can be performed. 



Nevertheless, a qualitative description of such a com- 
plex process can be constructed using simple ideas. Go- 
ing through the PDG listings, one reaUses that essentially 
the whole D decay width can be explained by simple tree- 
level valence quark diagrams, such as the one shown in 
Fig. [T] connected to the well known qq resonances from the 
Constituent Quark Model. Considering, for instance, three- 
body decays proceeding through intermediate states having 
spin-1 and 2 resonances, one concludes that the regular qq 
mesons correspond to the entire decay rate. No 'exotic' 
states have been observed in D decays, which act as a qq 
filter. Assuming the tree-level diagrams to be dominant, the 
available 'final state' quarks, that is, those resulting from 
the weak c decay, determine not only which qq resonances 
can occur, but also their relative rates. A nice example is 
the 'weak vector/axial-vector dominance', a close analogy 
to the vector dominance in electrodynamics. Due to the 
V-A nature of the c quark decay, in the diagram of Fig. [T] 
the virtual will couple preferentially to a vector or to 
an axial-veto particle, rather than to a pseudoscalar. The 
branching fraction for the decays D'^ — > (1260) and 

A'^/9+(770), for instance, are three to four times 
larger than that of £)° ^ if~7r+. 
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Figure 1: The T4^ -radiation diagram. In the case of the 
13+ i4r^7r+7r+ decay (q = d) the intermediate states 
are dominated by the K family (sd). 



On the analysis technique 

Essentially all studies of three-body hadronic D decays 
employ the same technique: the unbinned maximum like- 
lihood fit of the Dalitz plot, in which the decay matrix el- 
ement is represented by a coherent sum of phenomenolog- 
ical amplitudes [IJ. These amplitudes correspond to the 
possible intermediate states in the decay chain D Rh, 
R hh (h ^ K, tt). The amplitudes are grouped accord- 
ing to the orbital angular momentum L in the first step of 



the decay chain, 

L k 

The ampHtudes A]^ are weighted by constant complex 
coefficients c^, the series being truncated at L = 2. The 
set of complex coefficients is, in general, the fit output. 

In the case of a resonance with spin, the standard proce- 
dure is to define the resonant amplitude Ak as a product of 
a relativistic Breit-Wigner function, form factors (usually 
the Blatt-Weisskopf dumping factors |2|) for the D and the 
resonance decay vertexes and a function describing the an- 
gular distribution of the final state particles, accounting for 
the angular momentum conservation. The S-wave is the 
problematic issue. The way it is handled consists in the 
basic difference between the various Dalitz plot analyses. 
Here I will briefly describe the most common approaches 
to the S-wave. 

The Isobar Model 

In the so called Isobar Model the S-wave is usually as- 
sumed to be a sum of a constant nonresonant term and 
Breit-Wigner functions for the scalar resonances. The 
Breit-Wigner functions may or may not be multiplied by 
scalar form factors. In spite of conceptual problems, in 
most cases the Isobar Model provides a reasonably good, 
effective description of the data. 

In the beginning of this decade the pioneer work of 
the E79 1 Collaboration IH |5] showed evidence for two 
broad scalar resonances, identified to the a and the to k in 
the study of the — > 7r^7r"'"7r+ and D'^ — > if "tt+tt"'" 
decays, respectively. The E791 analysis used the Isobar 
Model, with one innovation: in addition to the constant 
complex coefficients c^, the masses and widths of the res- 
onances were also determined by the fit. 

The cr and the k, nowadays well established states, were 
soon after confirmed by other experiments, in different re- 
actions and with higher statistics ||6] |7] [H |9] • The values 
obtained by E791 for the Breit-Wigner masses and widths 
are inadequate for determining the <j and k poles. The 
merit of the £79 1 work was to demonstrate the existence of 
structures at low 7r+7r+ and K^-k^ mass with a resonant 
behaviour, that is, described only by an amplitude with a 
complex, energy-dependent phase. The Breit-Wigner was 
the simpler form of such an amplitude. 

In spite of yielding a good description of the data, the 
Isobar Model has a limited ability in disentangling individ- 
ual contributions from broad components in the S-wave. 
The case of the K^n^n^ decay is typical: the k 

and the nonresonant components are so highly correlated 
that the determination of the decay fractions become rather 
uncertain. 

This readily illustrated by the following exercise based 
on the result of the Isobar fit of K~-k'^tt^ Dalitz 

plot from FOCUS The set of coefficients from the 

FOCUS fit (Table II of ref. EOl is taken as die input model 



to simulate an ensemble with 2000 Dalitz plots. Each sim- 
ulated Dalitz plot had the same number (54K) of signal 
events as in the FOCUS data set. If there were no statistical 
fluctuations of the signal distribution, these 2000 samples 
would be identical. Each Dalitz plot was fitted with the 
same model. The resulting decay fractions were recorded. 
A scatter plot of the K7r+ and nonresonant decay fractions 
from the 2000 fits is displayed in Fig. |2l In the absence of 
correlations, the projection of the scatter plot onto each axis 
should look like a Gaussian centred at the value quoted in 
Table II of ifTOl ; the width should match the statistical error 
from the FOCUS data fit. But what we see is that the cor- 
relation between these two amplitudes is indeed very high, 
showing that the Isobar Model cannot provide a reliable 
distinction between the broad structures of the S-wave. 









□ □ □ □ = 












□ □□□□ 




□ □□□□ 




" □ □□□□ 




" ° □□□□ 




□□□□□ 












" □ □□□□ □ □ 




□ □□□□ 












□ =□□□□ = 




o □ □ □ □ 




□ o o □ □ o 






' . . . 1 . . 


. 1 ... 1 ... 1 ... 1 . . 1 . . . 



'20 22 24 26 28 30 32 34 36 38 40 



nonresonant decay fraction {%) 

Figure 2: The correlation between the K7r+ and the non- 
resonant amplitudes from the Dalitz plot analysis of the 
Z?+ K~iT^iT^ using the isobar model. See text for 
details. 



The K-matrix approach 

An alternative approach to the S-wave is the K-matrix 
formalism, applied to DaUtz plot analyses of D decays by 
the FOCUS collaboration ifTol lTll. This approach involves 
a very sophisticated machinery, but is based on a unrealis- 
tic and somewhat naive assumption: in the three-body fi- 
nal state, the tt'^tt^/K^tt^ pair forms an isolated system, 
which evolves as if the third body was not there. Three- 
body interaction is, therefore, ignored, as shown schemati- 
cally in Fig. [3] 

Given that a no rigorous treatment of a three-body fi- 
nal state strong interaction exists, this assumption greatly 
simplifies the problem. Only under this assumption one 
can invoke arguments such as two-body unitarity. In the 
absence of a full three-body final state interaction (FSI), 
the dynamics of the 7r^7r+7r+ and K^tt^tt^ final states 
becomes entirely determined by the two-body tt+tt^ and 
K^TT^ interactions, respectively. That allows one to con- 



strain the D decay amplitude by data from different reac- 
tions. The phase of the 7r+7r^/_ftr^7r+ ampHtude should, 
therefore, match that of the tt^tt^/K^tt^ scattering not 
only for the S-wave, but also for all other waves. That is 
the essence of Watson's theorem. 

There is no experimental evidence supporting this ap- 
proximation. As we will see, the S- and P-wave 
phases from D decays are rather different from that of 
T:+TT-/K-7t+ scattering lfm[T2l[T3l . 
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Figure 3: Schematic diagram of a three-body decay of a D 
meson in which the resonant system P1P2 does not interact 
with the third particle. This is the underlying picture in the 
K-matrix approach. According to it final state interaction 
occurs only between Pi and P2 ■ 

In the K-matrix approach the n^n^ IK^tt^ S-wave 
phase is fixed. The K-matrix approach is, hence, not in- 
tended to provide new data on the 7r+7r^/A'^7r+ S-wave 
phase. The decay amplitude is defined as a product of the 
K-matrix and a production amplitude. The fit parameters 
are contained in the production amplitude. That includes 
an adjustable energy dependent phase. In other words, the 
observed phase from D decays is modeled by the sum of 
the known phase from scattering and an unknown phase 
to be determined by the fit. The 'production phase' can, 
therefore, account for any eventual differences between the 
S-wave phase form tt^tt^ /K^ir^ scattering and from D 
decays. With such a freedom, the K-matrix approach pro- 
vides, in general, fits with acceptable quality. 

The MIPWA method 

The Model Independent Partial Wave Analysis technique 
was developed by the E791 Collaboration 1 11 1 , implement- 
ing an idea put forward by W. Dunwoodie. As in the iso- 
bar and K-matrix approaches, the decay matrix element is 
written as a sum of partial waves, truncated at the D-wave 
(which is already a very small contribution). No assump- 
tion is made on the nature of the S-wave, which is rep- 
resented by a generic complex function to be determined 
directly from data, 

Ao(s) = ao(s)e*^«(^). (1) 

The ■n^TT^IK^TT^ mass spectrum is divided into n slices 
in > 20, in general). For each slice two real numbers are 
fitted, so that at the k-th sUce Ao{s — Sk) = flgC^'^o. An 
interpolation is used to define the value of the S-wave in 
any point between Sk < s < Sk+i- The set of {a§, (/)§}, 
together with the coefficients are the fit parameters. 



In this sense, the MIPWA method is the most exempt 
approach. The only assumption is common to all other 
analyses, that is, the P- and D-waves are well described 
by a sum of Breit-Wigner amplitudes. The are some short- 
cuts, though. First, one should handle a large number of 
fit parameters (the S-wave alone has 2n free parameters), 
which introduces some technical difficulties. Moreover, 
the MIPWA S-wave relies on a precise representation of 
the other waves. If something is wrong with the P- and D- 
wave parametrisation, their content would "leak" into the 
S-wave. But the crucial problem is that the MIPWA S- 
wave is an inclusive measurement, since the tt+tt^/A' ^7r+ 
system is embedded in a three-body strongly interacting 
final state, as illustrated in Fig 2] Extracting the pure 
7r+7r~/Ar~7r+ amplitude is not a trivial task: one needs to 
deconvolute the desked phase from the ones introduced by 
three-body FSI lfT4l and, perhaps, from the production am- 
plitude. 




Figure 4: A more realistic schematic diagram of a Z) de- 
cay. Final state interactions may occur between all decay 
particles, including three-body interactions. 



The Kn amplitude - the k problem 

The existence of a broad scalar resonance at low K~Tr~^ 
mass was first reported by the £79 1 Collaboration |5| from 
the Dalitz plot analysis of the _D+ K^h^t:^ . This state 
was identified with the k(800) meson. Shortly after the 
same state was observed by several other experiments ||9] 

noiini. 

The the nature of the k(800) meson - an / = 1/2 state? 
- has been the subject of a long-standing debate. While 
there is now plenty of evidence for the neutral state, results 
for the charged partner are still scarce and conflicting lfT6l 
El- 

The k(800) pole position has been determined recently 
using LASS data 1 18] and a Roy-Steiner representation of 
K^n^ scattering amplitude 1 19|. Note, however, that there 
is no data on K^n^ K^tt^ bellow 825 MeV/c^, where 
LASS data starts. The crucial issues are, therefore: a) to 
fill the existing gap between the K^tt'^ threshold and 825 
MeV/c^; b) to find the chai-ged k(800) state. 

Searching for the charged k is a much harder task than 
for the neutral partner. There is no 'golden mode' in which 
the contribution of the {Ktt)^ in S-wave is largely dom- 
inant. The available data sets have still limited statistics. 



This will be illustrated by two studies, with a somewhat 
surprising results. 
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from BaBar 

The K^K^n^ decay was studied by BaBar 

Iflel . The Babar sample has 11000 events with 98% pu- 
rity. The K^K^tt'^ is a Cabibbo suppressed decay, 
with dominant tree-level amplitudes (external and inter- 
nal ly-radiation). The dominant contributions should come 
from the vector modes K*^K~ and 07r°, but we also ex- 
pect sizable decay fractions for the modes K*~ and 
{Kii)^K^. The i^TrS -wave can be analysed from thresh- 
old up to 1.4 GeV/c^ using the A'^tt^ system. 

The Dalitz plot of the — > K^K'^ir'^ decay is shown 
in Fig. |5j We see clearly the bands corresponding to the 
K*+K~, K*^K+ and cjiir^ modes (the KK axis runs 
along the top-right to bottom-left diagonal). The Dalitz plot 
projections are shown in Fig. |6l confirming the expectation 
of a larger K*^K^ contribution compared to K*^ . 
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Figure 5: The Dalitz plot of the D° K-K+tt° decay. 
The narrow band on the top right part of the plot corre- 
sponds to the 07r° mode. The horizontal structure is due 
to the dominant mode, the K*^K^ . The mode K A'+ 
appears as a vertical band. In all cases we see the node due 
to the angular distribution typical of a spin-1 resonance. 

The Dalitz plot of Fig. 5 was fit with three different mod- 
els for the K^tt'^ S-wave. In the first model, the S-wave 
was represented by the LASS 7=1/2 amplitude. 



Aq(s) = ^sin(5(s)e"^('') 
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S{s) = cot 
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In the above equation s is the Kir mass squared, a and 
b are real constants and p is the kaon momentum in the 
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Figure 6: Projection of the D'^ K~ K'^ir^ Dalitz plot of 
Fig. |5] The projection on the K^tt^ axis (left plot) show 
the peak of the K* (892)+, much more prominent than the 
one in the K 7r° projection (right plot). In both plots the 
structure between 1-2 GeV^/c^ is the projection of the 07r° 
mode. 

Ktt rest frame. The parameters Mq and Fq refer to the 
Aro(1430) resonance. 

In the second model BaBar used the E791 MIPWA 
K~'K^ amplitude fTTl. In the third model a coherent sum 
of a uniform nonresonant term and Breit-Wigner ampli- 
tudes for the Tfo"^ (1430) A:=f and k(800)± AT^ modes - the 
isobar model. 

The isobar model yielded the smaller fit probability. The 
best fit was obtained with the LASS 7=1/2 S-wave am- 
plitude, although a good fit was also achieved with E791 
MIPWA S-wave. The latter describes the data well, except 
in the region near threshold. 

BaBar data is well described by two models, both using 
the LASS 7=1/2 amplitude for the K^ir^ S-wave. The first 
model has nine amplitudes, whereas the second has only 
six. The basic difference between the two models are the 
74r*(1410)^7C+ modes, present in the first model and ab- 
sent in the second one. The decay fractions from the two 
fits are listed in Table [T] 

The failure of the isobar model in describing the data 
cannot be taken as an argument against the k. Recall that 
the K pole was found in LASS data. One would expect the 
LASS amplitude to be well suited for situations where the 
Ktt system is isolated from the rest of the final state, as 
in semileptonic decays. The fact that the LASS amplitude 
yielded the best fit is a bit surprising. 

We should analyse this result with some care, though. 
An inspection of the BaBar fit fractions in Table I shows 
that the interpretation is not straightforward. The decay 
fractions from model II, where the tiny contributions of 
the K*[l'ilO)^K^ modes were removed, add up to 165%! 
The fraction of the 7C*(1410) resonance in D decays is al- 
ways marginal, as in the present case. However, when this 
small component is removed from the fit the P-wave re- 
mains unaltered, but the TT+tt" S-wave contribution jumps 
from 16.3% to 71.1%. The TC^tt" S-wave component, con- 
sistent with zero in model I, becomes a lOu effect in model 
II. There is an obvious interplay between the 7C*(1410)7C 
and the TC^tt" S-wave. With more data the TiT+Tr" S-wave 
amplitude could be extracted with the MIPWA technique. 
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K*{892)+K- 


45.2 ± 0.9 


44.4 ± 


0.9 


K*{141Q)+K^ 


3.7± 1.5 








16.3 ±0.1 


71.1 ± 


4.2 




19.3 ± 0.7 


19.4 ± 


0.7 


/o(980)7r" 


6.7 ± 1.8 


10.5 ± 


1.4 


K*{892)-K+ 


16.0 ±0.9 


15.9 ± 


0.9 


K*iU10)-K+ 


2.7± 1.5 








2.7± 1.5 


3.9 ± 


1 



Table 1: Decay fractions, in %, from the BaBar Z?" 
K^K+TT° Dalitz plot fit using the LASS 7=1/2 amplitude 
for the X='=7r° S-wave . As expected, the K*{892)+K^ 
is the dominant mode. Removing the small K*{14:10)K 
component causes a drastic change in the K'^tt'^ S-wave 
contribution, bringing the sum of the decay fractions to 
over 165%. 

r~ ^IC^TT^Ur from Belle 

Semileptonic decays like D Kiilv and r — > Kttv 
are very interesting alternatives, since the Kit system is 
free from final state strong interaction. We should expect 
Watson's theorem to hold, or, in other words, that the Kit 
S-wave phase matches that from LASS. There are some 
problems with semileptonic decays, though. In these de- 
cays the P-wave corresponds to over 90% of the decay rate. 
Very large samples are required in order to have a reason- 
able statistics for the S-wave. 

The ^ K TT^Vr decay was studied by Belle |17|. 
The sample was selected from events of the type e+e~ — ^ 
T+r^, with T+ l^VrVi and Ksir^Vr- The sig- 

nature was a lepton recoiling against a pair of pions of op- 
posite charge. The selected sample has 53K signal events. 

The Ks7T~ mass spectrum is shown in Fig. |2l Super- 
imposed (histogram in red) we see the K (892)- con- 
tribution. There is an excess of data events over the 
K (892)~ contribution both at the lower and the higher 
part of the KsTt~ spectrum. The spectrum of Fig. [7] was 
fitted with different models. To the dominant K (892)- 
two other amplitudes were added: the k(800)~ plus one 
Ktt resonance with higher mass — either the K (1410)-, 
i?o(1430)- or the if* (1680)-. The LASS amplitude (eq. 
IJlandO was also tried. 

The result of the fit of the KsTt^ spectrum was also sur- 
prising. Contrarily to what one would expect, the model 
with the LASS 7=1/2 amplitude fails to reproduce the 
KgTT^ line shape (C.L.=10-*). The best description of the 
data was achieved by adding to the K (892)- a pure scalar 
component, that is, the model with the k(800)~ plus the 
7^0(1430)- resonance. 

The missing neutrinos introduce additional limitations. 
The full event reconstruction becomes very difficult. One 
has to handle a relatively high background, at the 20% level 
in this analysis. The most serious consequence is that, since 
the position of both the primary and secondary vertexes are 
not determined, no angular analysis can be performed. 




Figure 7: The TC^tt spectrum from the r K ir Vr 
decay, from Belle. The solid histogram is the contribution 

of the 7?* (892)-. 

In the case of the 7?+ K^n^ fi^v decay from FO- 
CUS II20I the angular distribution was a crucial piece of 
information. In this decays there is a 7% contribution from 
the 7f-7r+ S-wave. The S-wave component interferes with 
the K (892)~, causing an asymmetry in the helicity angle 
distribution (the helicity angle is defined as the angle be- 
tween the kaon momentum and the line of flight of the 7)+, 
measured in the K^tt^ rest frame). The K^n^ line shape 
from 7)+ K^TT^ ^'^v could be fitted equally well with 
different S-wave models, but each model has a different 
interference pattern with the P-wave, distorting the helicity 
angle distribution in a different way. The distribution of the 
helicity angle could then be used to discriminate between 
the different S-wave models. 

A strong case for the k(800)- would be made from the 
^ K 'k~Vt decay if the angular analysis was per- 
formed. Unfortunately no information on the angular dis- 
tribution is available. 

K-TI+TI+ from FOCUS 

Let's now turn to a situation where the Ktt S-wave is 
largely dominant. The issue here is the S-wave phase 
near threshold, which can only be addressed by heavy fla- 
vor decays. The 7)+ K^tt^tt^ is a 'golden mode': 
large branching fraction, easy to be reconstructed, very low 
background and with an S-wave contribution amounting to 
approximately 80% of the total decay rate. 

This decay was studied in great detail by FOCUS ifTOl 
[T2I . A sample with 54K signal events and 98.5% purity 
was analysed with the MIPWA technique. The 7)+ ^ 
7C~7r+7r+ Dalitz plot is shown in Fig. |8] Since there are 
two identical pions, the Dalitz plot is symmetric. The blue 
lines indicates the 7ir*(892) mass squared. We see clearly 
the effect of the angular distribution splitting the 7ir*(892) 
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Figure 8: The D+ K-tt+tt+ Dalitz plot from FOCUS. 
Due to the identical pions, the plot is symmetric with re- 
spect to the diagonal (TO^(7r7r)) axis. 



band into two lobes. A striking feature is the displacement 
of the two lobes with respect to the nominal K* (892) mass. 
This is readily explained by the interference between the 
K (892)7r+ mode and the S-wave. This interference al- 
lows one to measure the S-wave phase. 
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Figure 9: The K 7r+ MIPWA S-wave (circles with error 
bars) from FOCUS D+ K-ti+tt+ decay. The LASS 
1=1/2 S-wave phase ((5i/2) is shown as full stars. The Ktt 
amplitude is elastic up to 1.45 GeV/c^, indicated by the 
vertical Une. 



ii'2(1430)7r+ is also present. 

The FOCUS MIPWA S-wave phase is shown in Fig. |9] 
The dotted line indicates the Kij' threshold, up to which 
the Ktt scattering amplitude is elastic. The circles with 
error bars are the FOCUS result, whereas the black stars 
are the LASS 7=1/2 S-wave phase. The elastic region is 
highlighted in Fig. [TOl which shows the FOCUS S-wave 
phase, shifted by 80°, together with the 7=1/2 and 7=3/2 
S-wave phases from LASS. 

All resonances are in the 7=1/2 component, while the 
7=3/2 amplitude is purely nonresonant. Since in D de- 
cays the nonresonant contribution is usually very small, one 
would expect the FOCUS phase to be similar to the LASS 
7=1/2 phase. But we see that this is not the case. As a 
matter of fact, no combination of the two LASS isospin 
components can reproduce the S-wave phase from FO- 
CUS. The S-wave phases from 73+ 7C^7r+7r+ and from 



K- 



K 7r+ are indeed very different. An additional 



energy dependent phase must be added to the LASS phase 
in order to match the FOCUS result. 

Why are the two phases so different? Where this ad- 
ditional energy dependent phase comes from? There are 
two possible origins: the decay amplitude and three-body 
final state interaction. It is rather suggestive that the differ- 
ence between LASS and FOCUS phases increases as one 
approaches the Kn threshold. As the Ktt mass decreases, 
the momentum of the third particle increases, and its in- 
teraction with the Ktt system becomes more intense. That 
would explain why the E791 MIPWA S-wave amplitude 
does not yield a good fit to the BaBar K^n^ 
Dalitz plot at low Ktt masses: at the TC+tt*' threshold, the 
third particle (a K^) has smaller momentum, so the phase 
introduced by the three-body FSI would be slightly differ- 
ent than that from the 7)+ K~tt^tt^. 
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The S-wave component corresponds to 80.2% of the 
decay rate. The remaining part is mostly due to the P- 
wave, which is well described by the K (892)7r+ and the 
!K'*(1680)7r+ modes (the contribution of the !^*(1410)7r+ 
mode is consistent with zero). A small fraction of 



Figure 10: The FOCUS K-7r+ MIPWA S-wave shifted by 
80°. The lines are the LASS isospin 1/2 and 3/2 ampli- 
tudes. Only the elastic region is shown. 

The basic problem with the MIPWA S-wave is, there- 



fore, how to interpret the data. The pure K-k amplitude is 
there, covering the entire elastic range from threshold. Un- 
fortunately it is not directly accessible. We need to learn 
how to get it. 

The TTTT amplitude - the /o(1370) problem 

The TTTT S-wave in the region 1.2-1.5 GeV/c^ is still prob- 
lematic. There are two states in this region, namely the 
/o(1370) and /o(1500). The /o(1500) is a well estabUshed 
resonance, observed clearly in pp and pp data, and also in 
J /if) decays. Its mass is (1.505±6) GeV/c^ and its width is 
(109±7) GeV/c'^, having also well measured couplings to 
TTTT, Att, KK and rjri UJ. 

The /o(1370), by its turn, remains very controversial. Its 
mass ranges from 1.2 to 1.5 GeV/c^, while the width lies 
between 200 and 500 MeV/c^ H]. The BES Collaboration 
observed an excess of events in 1.2-1.5 GeV/c^ region in 
the decay J /if) (/)7r+7r~ II2TII . which was interpreted as 
a dominant /o(1370) component interfering with a small 
contribution of the /o(1500). The /o(1370) was repre- 
sented by a Breit-Wigner and the values obtained for the 
mass and width were (1.350±50) and (0.265±40) GeV/c^, 
respectively. No evidence of the /o(1370) was found in 
J/^ (t>K+K- and J/t/i -f-K+n- /-iK+ R- . The 
ratio of partial widths obtained by BES is consistent with 
zero: Vkk/^^-,.= (0.08±0.08). 

The region around 1.5 GeV is very interesting: that's 
where the ground state of the scalar glueball is expected 
to be. It is necessary, therefore, to measure not only the 
/o(1370) mass and width, but also the couplings to other 
channels, for this would provide insight to its nature. One 
possible scenario includes also the /o(1710). The three ob- 
served states would be mixtures of two qq and the 0++ gg 
states [22J. 

The information given by heavy flavor decays is partic- 
ularly useful in this respect. The states that are observed 
with a large decay fraction in D and B decays are very 
likely to have a dominant qq component. 

Dt 7r-7r+7r+ from FOCUS and E791 

The D'^ 7r^7r+7r+ is a 'golden mode' for studies of 
the TTTT system in S-wave. This is a Cabibbo suppressed 
mode with no strange quarks in the final state. Resonances 
that couple both to KK and to tttt, like the /o(980), are 
expected to play a dominant role. 

This decay was studied by E791 |3 1, and, more recently, 
by FOCUS |71. In Fig. [lUwe see the D+ tt-tt+tt+ 
Dalitz plot from FOCUS. Two features call the atten- 
tion immediately: the narrow bands corresponding to the 
/o(980)7r+ mode and the concentration of events between 
1.5-2.2 GeW^/c'^. This concentration is partially due to a 
scalar state with high mass, which will referred to as the 

MX). 

A Dalitz plot analysis was performed both with the K- 
matrix and with the isobar model (but only the result of K- 
matrix fit were published). In addition to the /o(980), the 



3.5 




a F , , , , I , , , , I , , , , I , , , , "r, : , , I , , , , I , , , 
0.5 1 1.5 2 2.5 3 3.5 

Figure 11: The 13+ 7r^7r+7r+ Dalitz plot from 

FOCUS. The narrow bands at 1 GeV/c^ correspond to 
the /o(980)7r+ channel. The concentration of events at 
2 GeV/c^ are due to a scalar state identified with the 
/o(1500). 

S-wave isobar model included a scalar state, the fa{X), 
represented by a relativistic Breit-Wigner. The mass and 
width of the fo{X) were determined by the fit. The values 
obtained by FOCUS are (1.476±6) GeV/c^ for the mass 
and (0.119±18) GeV/c^ for the width. The data is well 
described by an S-wave with only two resonances. 

Monte Carlo simulations of the decay Df — > /o(X)7r+ 
were performed taking the values obtained by BES and FO- 
CUS for the mass and width of the foiX). For compari- 
son, simulations were performed assuming for the fo{X) 
the PDG values for the /o(15GO) and /o(1710). The sim- 
ulations are shown in Fig. [12] The conclusion is that the 
scalar state observed in FOCUS data is much closer to the 
/o(1500) than to the state observed by BES. 

I}+ ^ 7r~7r+7r+ from BaBar 

Recently BaBar reported study of the D+ 7r~7r+7r+ 
decay [13], from a high purity sample of about 13000 
events. The 7r^7r+ S-wave was measured using the 
MIPWA technique. The result is displayed in Fig. [13] In 
the left plot we see the S-wave magnitude as a function of 
the TT+TT^mass. There are two peaks, a narrow one at the 
/o(980) mass and another at 1.4-1.5 GeV/c^, which is rel- 
atively narrow. In the plot on the right we see the S-wave 
phase, also as a function of the 7r+7r^mass. There is a rapid 
variation of the phase as one crosses the /o(980) mass, as 
expected for a typical resonance behaviour. The phase con- 
tinues to grow and between 1.4-1.5 GeV/c^ another rapid 
variation can be observed, indicating the presence of an- 
other resonance. The magnitude and phase of the S-wave 
from the FOCUS analysis is superimposed to the BaBar 
result. The agreement between FOCUS and BaBar is very 




fo(1370) 
(BBS) 



3.5^ 

N 3 ; 
2.5; 

2; 

1-5 Jiiii 

1 M 
0.5 ; 





fo(1500) 
(PDG) 




FOCUS j 


(C) J 

- 


^ E79I 1 
- i 












m(j:+ 71') (GeV/c-) 



m(7C* Til (GeV/c^) 



fo(1475) 
(FOCUS) 




Figure 12: Monte Carlo simulation of the Dalitz plot of 
the Df /o(X)7r+ decay, with different hypothesis for 
the foiX). Clockwise from top left: /o(1370) (BBS), 
/o(1500) (PDG), /o(1710) (PDG) and the state found in 
FOCUS data. 



good. 

Although the conclusion drawn from Fig. [T2]is that the 
foiX) is consistent with the /o(1500), the values of the 
mass and width obtained from the Df 7r~7r+7r+ Dalitz 
plot fit are not quite the same as the PDG values for the 
/o(1500). We should keep in mind that the Breit-Wigner 
which was used is only an approximate representation for 
this state. In FOCUS analysis the tttt mode was assumed 
to account for all the /o(1500) decay rate. The total decay 
width should be T{s) = r^^(s) + T4^{s) + Tkk{s) + 
r,pj(s). Moreover, since the Z?+ mass is not too high, the 
/o(1500) peak lies out of the Dalitz plot boundary. 



from Belle 



Charmless three-body B decays are a very promising 
tool for light quark spectroscopy. As in the case of charm 
decays, charmless B decay have a rich resonant structure. 
The phase space of B decays is much larger than that of 
D decays, so resonances are fully contained in the Dalitz 
plot. However, since the branching fractions are typically 
between 10~^-10~^, the statistics is still limited. This will 
no longer be an issue when the LHCb data becomes avail- 
able. 

Two decay modes are particularly interesting for the 
/o(1370) problem: B+ -> K+tt+tt- and B° 
The dominant mechanisms in the B^ 



A TT^ 

K TT+TT^ decay are assumed to be the penguin diagram of 



Figure 13 



Results from the BaBar MIPWA fit of the 
+ Dalitz plot, showing the magnitude (left 
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plot) and phase (right plot) of the it~ 
tion of the tt^tt^ 
posed to the BaBar result show the magnitude and phase of 
the FOCUS and E791 S-wave form the isobar fit. 



Fig. [14] and a tree-level Cabibbo suppressed diagram (ex- 
ternal W-radiation). The diagrams for 5+ A'+tt+tt" 
are obtained replacing the d by an u quark. One ex- 
pects, therefore, the same intermediate states, except for 
the charge of the Ktt resonances, with similar decay frac- 
tions. 




Figure 14: The dominant amplitude for the B^ 
TT~ decay. The diagram for the _B+ K^n^ 



is obtained replacing the d quarks by u. 

These two modes were studied by Belle 12311241 . In Figs. 
[15] and [16] we see the projection of the i?+ — > K^tt^tt~ 
Dalitz plot onto the K^tt~ and tt+tt^ axes, respectively. 
The K^TT^ projection exhibits two prominent structures, 
the narrower corresponding to the A'*(892)7r+ decay and a 
broader corresponding to the A'q (1430)7r. 

In FigfTSl three peaks are clearly visible, corresponding 
to the p(770), to the /o(980). The third peak is well de- 
scribed by a model with only one high mass scalar state at 
-1.4-1.5 GeV/c2. 

The decay fractions from the Belle Dalitz plot fits 
are presented in Table |2] The empirical parametrisation 
aie^^ e""*^'" +026*^ g-as„„ was used for the nonresonant 
amplitude. The nonresonant contribution is dominated by 
the Ktt component in both B+ and decays. The very 
high fraction of the ATq (1430)7r fraction is a bit surprising. 
It is hard to believe that the ArQ(1430)7r decay fraction is 
five to six times larger than that of the Ar*(892)7r. Here 
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Figure 15: The tt+tt projection of the £?+ K^tt^tt 
DaUtz plot from Belle. 
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Figure 16: The K^tt projection of the _B+ K^tt^tt 
Dahtz plot from Belle. 

we see again the interplay between broad structures in the 
S-wave, resulting in abnormally large decay fractions. 

However, when we look to the tttt component, every- 
thing seems under control. The peak at 1.4-1.5 GeV/c^ 
in the tt+tt^ projection is well described by a sin- 
gle scalar resonance, modeled by a Breit-Wigner am- 
plitude. The Breit-Wigner parameters of this resonance 
were obtained from data: Ma =(1.449±0.013) GeV/c^, 
To =(0.126±0.025) GeV/c^. These values are in good 
agreement with the ones from FOCUS. 

It seems that in heavy flavor decays only one scalar state 
is observed in the tttt channel. This state is not consistent 
with the /o(1370). It is similar to the /o(1500), although 
the masses and widths obtained from Dalitz plot analysis of 
D and B decays are not quite the same. The difference is 
not large and could be attributed to the way the parameters 
were determined by FOCUS and Belle. In any case, here- 
after we will refer to this scalar resonance as the /o(1475). 



mode 


B+ K+TT+n- 




K*{892)tt 


13.0±1.0 


11.8±1.7 


ii:o*(1430)7r 


65.5±4.5 


64.8±7.8 


p{770)K 


7.9±1.0 


12.9±2.0 


fc)(980)K 


17.0±3.6 


16.0±4.2 


foiX)K 


4.1±0.9 


3.7±2.4 


nonresonant 


34.0±2.7 


41.9±5.5 



Table 2: Decay fractions from the Belle B Kirir Dalitz 
plot fits. The fractions for both modes are in good agree- 
ment, as expected from an isospin symmetry argument. 



D+ ^ K+K-TT+ from CLEO-c 

One interesting aspect is the large decay fraction of the 
/o(1475) observed in the £)+ — > 7r~7r+7r+. Assuming the 
main decay mechanism to be the VF-radiation (Fig. 1), the 
large fraction of the /o(1475) in 7r^7r+7r+ may 

be interpreted as an indication of a strong ss component in 
wave function of this state. In this case, a large contribution 
would also be expected in _D+ — > K'^K~'k^ . 

The Z?+ K^K^TT^ decay was studied by CLEO-c 
125). The CLEO-c sample has 14K events with very small 
background. The Dalitz plot is shown in Fig. [17] This is 
a very tough analysis. The lower part of the KK spectrum 
is populated by the the /o(980) and the ao(980). There 
is a strong interference between these amplitudes, so it is 
very difficult to separate individual contributions. We see 
also the bands corresponding to the K*{892)K^ mode. In 
the region mj^j^ ^2.Q-2.2 GeV/c^, however, there is no 
indication of a resonance. 

This is very intriguing. Apparently the /o(1475) does 
not decay to KK. This is in agreement with the /o(1500) 
pai-tial width TKK/Tioi = 0.086 ± 0.010, and may be 
considered as an additional evidence for the identification 
of the two states. But the mechanism that leads to a large 
/o(1475) decay fraction in _D+ 7r^7r+7r+ remains to be 
understood. 

Conclusions 

The physics of the scalar mesons has still many inter- 
esting open problems, which are related to the strong dy- 
namics at low energy. The main problem in this field is 
to identify the qq mesons of the Constituent Quark Model 
scalar nonet(s). There are currently more candidates than 
slots, although some states remain controversial. From the 
experimental point of view, it is not trivial to detect broad, 
structureless overlapping states squeezed in a limited phase 
space, a situation which is well illustrated by the Colom- 
bian painter Fernando Botero, in Fig [18] In addition to 
the regular qq states, there is a number of other configura- 
tions allowed by QCD, like molecules, hybrids, glueballs, 
tetraquarks, sharing the same J^*^ ~ 0++ quantum num- 
bers. These 'exotic' configurations have not been clearly 
identified yet, but they may be mixed with the regular qq 
mesons. 
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Figure 17: The D+ K+K-'k+ DaUtz plot from CLEO- 

ciEa. 



It is, therefore, unlikely that the understanding of the 
nature of scalar particles could be achieved without com- 
bining data from different types of reactions. Heavy fla- 
vor decays have been explored as an alternative window 
to some fundamental issues, like the nature of the tttt and 
Kit spectrum near threshold. In this work the k(800) and 
the /o(1370) were discussed from the point of view of 
hadronic three-body decays of D and B mesons and also 
from decays of t lepton. 

The neutral k(800) is well established. Its pole posi- 
tion was determined, in spite of the lack of data bellow 
825 Mev/c^ in the K~ti^ spectrum. This gap can be filled 
by data from _D+ K^tt^tt^ using the MIPWA tech- 
nique, but we need to understand what exactly is being 
measured, how to account for three-body final state inter- 
actions, whether or not the decay dynamics introduces an 
energy dependent phase. 

If the k(800) is a 7=1/2 state, then its charged part- 
ner must exist. The search for the k(800)^ is a hard 
task, though. The cleanest environment is provided by the 
semileptonic decays of D mesons, such as 73 ^ K-Kjiv. 
However, very large samples are required, since the S-wave 
is only a small component. Decays of t leptons would 
be an interesting alternative, but the missing neutrinos is 
a serious obstacle. Another alternative is to measure the 
S-wave from decays like D° K~K+tt° (LHCb). The 
problem here is twofold. In addition to the same difficulties 
as in the 73+ 7i ^7r+7r+, in hadron machines it is harder 
to select a clean sample of modes with neutral pions. At 
this point, the nature of the k remains an open question. 

The situation concerning the /o(1370) is also a bit ob- 
scure. The existence of this state might not even be taken 



for granted. If it had a large qq component, it should 
have been unambiguously observed in charm decays. The 
scalar state that is present in D and B decays is, in- 
stead, similar to the /o(1500). The MIPWA analysis of 
the 73+ 7r^7r+7r+ show a resonant behaviour between 
1.4-1.5 GeV/c^, in agreement with FOCUS and E791 find- 
ings. The Breit-Wigner parameters of this scalar state are 
not quite compatible with the PDG values for the /o (1500). 
One important aspect that favours the identification of this 
scalar meson with the /o(1500) is that it is not seen in 
73+ 7C^7f+7r+. With more data and a refined anal- 
ysis technique the pole position of the /o(1475) could be 
determined. 

Very soon we will enter the era of 'infinite statistics' in 
essentially all interesting decay modes. The expected num- 
bers from LHCb are really impressive. But we must ac- 
knowledge that today we are already limited by systemat- 
ics. The experimentalists are not ready to explore the full 
potential of the coming data. The models currently used 
to parametrise the signal distributions, like the Dalitz plot, 
are inadequate. New analysis tools with a better theoretical 
foundation are urgently necessary. 




Figure 18: Fernando Botero: Dancing in Colombia 
(1983.251). In Heilbrunn Timeline of Art History. 
New York: The Metropohtan Museum of Art, 2000. 
[http : //w w w. me tmuseum .org/toah/ho/ 11/ sa/ho- 1 9 8 3 . 25 1 .htm 
(October 2006). 
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